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Abstract 

This paper exploits the connection between the quantum many-particle density of states and 
the partitioning of an integer in number theory. For N bosons in a one dimensional harmonic 
oscillator potential, it is well known that the asymptotic {N — > cxd) density of states is identical to 
the Hardy-Ramanujan formula for the partitions p{n), of a number n into a sum of integers. We 
show that the same statistical mechanics technique for the density of states of bosons in a power- 
law spectrum yields the partitioning formula for p^{n), the latter being the number of partitions 
of n into a sum of ■s*'' powers of a set of integers. By making an appropriate modification of the 
statistical technique, we are also able to obtain d^{n) for distinct partitions. We find that the 
distinct square partitions d'^{n) show pronounced oscillations as a function of n about the smooth 
curve derived by us. The origin of these oscillations from the quantum point of view is discussed. 
After deriving the Erdos-Lehner formula for restricted partitions for the s = 1 case by our method, 
we generalize it to obtain a new formula for distinct restricted partitions. 
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INTRODUCTION 



The N-particle density of states of a self-bound or trapped system has attracted the 

attention of physicists for a long time. We have in mind the work done in nuclear 111 and 

□ 1 

particle physics as well as in connection to black-hole entropy 3| in recent times. In 

nuclear physics, one is generally interested in self-bound fermions at an excitation energy 
E that is large compared to the average single-particle level spacing, but is small compared 
to the fermi energy of the nucleus. In this energy range, the density of states is given by 
the highly successful Bethe formula ^ that grows as exp{ay/E), and is insensitive to the 
details of the single-particle spectrum. The constant a in the exponent is proportional to 
the single particle density of states at the fermi energy Ep. In hadronic physics l^, the 
many-particle density of states grows exponentially with E, and leads to the concept of a 
limiting temperature. The same behavior is found to hold for a bosonic system like gluons 
in a bag Q| . 

It is well known that for ideal bosons in a one- dimensional harmonic trap, the asymptotic 
( ^ oo ) density of states is the same as the number of ways of partitioning an integer n 
into a sum of other integers, and is given by the famous Hardy- Ramanuj an formula p. It 
also grows exponentially as VE, the same as the Bethe formula when E is identified with 
n. Grossmann and Holthaus lol have studied this system, and have used more advanced 



results from the theory of partitions J71| to calculate the microcanonical number fluctuation 
from the ground state of the system as a function of temperature. Combinatorial methods 
have also been used to compute the thermodynamic functions for similar systems ^. In 
this paper we use the N-particle quantum density of states (that may be derived using the 
methods of statistical mechanics) to obtain some novel results on the partitioning of an 
integer into a sum of squares, or a sum of cubes, etc. To the best of our knowledge, some 
of the results pertaining to the partitions of an integer to a sum of distinct powers are new, 
and will be pointed out as they appear in the text. For the harmonic spectrum, we are also 
able to obtain the leading order finite-N ( Erdos-Lehner [3| ) correction to the asymptotic 
Hardy- Ramanuj an formula using our method, and then get the corresponding (new) result 
for distinct restricted partitioning. 

In section II of this paper, we consider ideal bosons with a single particle spectrum given 
by a sequence of numbers generated by m^{m = 1, 2, 3, ■ ■ ■) for a given integer s > 1. We 
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first derive the asymptotic many-particle density of states for this system using the canonical 
ensemble and in the saddle-point approximation, and show that it grows exponentially as 
the (s + l)th-root of the excitation energy. Specifically, in the physically relevant case of a 
square well, this result implies that the density of states grows exponentially as the cube 
root of energy. Our general expression for the asymptotic density of states agrees with the 
Hardy-Ramanujan formula 9] for p'^{n), the number of ways an integer n may be expressed 
as a sum of s*'' powers of integers. Throughout this paper, we drop the superscript s when 
s = 1. 

The statistical mechanics method can also be modified to obtain asymptotically the num- 
ber of distinct partitions d'^{n) of an integer n using the partition function of the fermionic 
particle spectrum (excluding the hole distribution, to be explained later). This analysis is 
presented in section III, where the smooth part of the asymptotic density of states (which 
reproduces the average behavior of the distinct partitions) is derived using the saddle-point 
approximation. Interestingly, for the s = 2 case where the integer n is expressed as a sum of 
distinct squares of integers, computations of the exact values of (P{n) reveal large fluctuation 
about the smooth average curve. These fluctuations wax and wane in a beating pattern. The 
ratio of the amplitude of the oscillations to the smooth part of (P{n) goes to zero as n ^ oo. 
From the quantum angle, these oscillations in the many-particle density of states have their 
origin in the fluctuation of the degeneracy of the many-particle density of states about the 
average value. This, in turn, is related to the oscillatory part of the single-particle density 
of states of a one- dimensional square-well potential, and the constraint brought about by 
the Pauli principle. 

In section IV, we discuss the corrections to the saddle-point approximation when the 
number of particles is finite. In the theory of partitions this is known as restricted 
partitions (because of the upper limit on the number of partitions), as opposed to the 
unrestricted partitions discussed in sections II and III. We restrict our derivations to the 
harmonic oscillator (s = 1) spectrum in this section, since the canonical partition function 
is exactly known even for finite N in this case. For the bosonic case, using this partition 
function, we derive exponentially small corrections for finite A^, a result that agrees with the 
Erdos-Lehner [3| asymptotic formula for s = 1. This method is then extended for finding 
the finite-N correction for distinct partitions, a result that to our knowledge is new. We 
conclude with a summary of the main results, as well as some limitations of the present 
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work. 



THE MANY-PARTICLE DENSITY OF STATES 

We first discuss the general statistical mechanical formulation for a A^— particle system. 
The canonical A^— particle partition function is given by 



Zn{I3) = E m^M-^Ef^) = r PN{E)exp{-f3E)dE , (1) 

i 

where f3 is the inverse temperature, E^^^ are the eigenenergies of the A^— particle system 
with degeneracies r/j, and Pn{E) = J2i Vi^{E — E\^^) is the A^— particle density of states. The 
density of states Pn{E) may therefore be expressed through the inverse Laplace transform 
of the canonical partition function 

1 rioo 

Pn{E) = — / exp{f3E)ZN{(3)d(3 . (2) 
zm J -too 

In general, it is not always possible to do this inversion analytically. Note that the single- 
particle density of states may be decomposed into an average (smooth) part, and oscillating 
components This, in turn, results in a smooth part 'p^[E)^ and an oscillating part 
6pn{E) [3] for the A^— particle density of states : 

p^iE)=p^iE) + 6p^iE) . (3) 

The smooth part Pi^{E) may be obtained by evaluating Eq. © using the saddle-point 
„.thod Q. Unlike the one-particle case, where the oscillating part n,ay be obtained n.ng 

the periodic orbits in a "trace formula" |3] , it remains a challenging task to find an expression 
for the oscillating part 5pn{E) 10]. In what follows, we shall use the saddle-point method 
to obtain the smooth asymptotic pj^{E), and identify it with the Hardy-Ramanujan formula 
for p^{n). 

Before doing this, we note that the canonical partition function Zjs[{(3) for a set of non- 
interacting particles with single-particle energies e^, occupancies {nj}, may also be written 

as 

Z^{I3) = exp(-/3Ef Q{N, E,) exp{-f3E,{n,}). (4) 

{Ui} 

In the above, -E'o^^ is the ground-state energy which we set to zero, and E^ is the excitation 
energy. The sum is over the allowed occupation numbers for particles such that E^ = J2i ^i^i- 



Note that for a given E^, the number of excited particles in the allowed configurations may 
vary from one to a maximum of N. We denote by ^(A^, E^) the total number of such distinct 
configurations allowed at an excitation energy E^. We set the lowest single-particle energy 
at zero in order that Eq'^^ = 0, and consider a single-particle spectrum e^n = m^- If now the 
excitation energy E^ takes only integral values n, then Q{N, E) is the same as the number of 
restricted partitions of n, denoted by p%{n), and asymptotically equivalent to the density of 
states 'pj^{E). Omitting the subscript A^, as in p^{n), will imply that we are taking N oo, 
corresponding to unrestricted partitioning. 

To perform the saddle-point integration of Eq. 0, note that the integrand may be written 
as exp[S{P)], where S{(3) is the entropy given by, 

S{(3) = (3E + \ogZN. (5) 

Expanding the entropy around the stationary point /Sq and retaining only up to the quadratic 
term in the expansion in Eq. Q yields the standard result 

. -£E(M. (6) 

where the prime denotes differentiation with respect to inverse temperature and 

We now proceed with a single particle spectrum given by e^, = m**, where the integer 
m > 1, and s > for a system of bosons. The energy is measured in dimensionless units. For 
example, when s = 1 the spectrum can be mapped on to the spectrum of a one dimensional 
oscillator where the energy is measured in units of hu. For s = 2, it is equivalent to setting 
energy unit as /2m, where is m is the particle mass in a one dimensional square well 
with unit length. These are the only two physically interesting cases. We however keep s 
arbitrary even though for s > 2 there are no quadratic hamiltonian systems. In particular 
s need not even be an integer except to allow a comparison between the number theoretic 
results for p%{n) and the density of states Pn{E) that we obtain here. We first obtain 
the asymptotic results for unrestricted partitioning by letting ^ oo and discuss the N- 
dependent correction later. The canonical partition function in this limit may be written 
as 

oo -j^ 



where we have used the power-law form for the single particle spectrum. By setting x = 
exp(— we see that the bosonic canonical partition function is nothing but the generating 
function for p'^(n) in number theory, the number of partitions of n into perfect s-th 
powers of a set of integers : 

OO OO -1 

zM = Y.p'inK=llTr—:^y (9) 

n=l n=l \ 

In the limit N oo, p^{n) is the same as Q{E) where the energy E is replaced by the integer 
n. In general the above form holds for all s in the limit of — oo, but is exact for finite N 
only for the oscillator (s = 1) system Using Eqs. (jSllHl), and the Euler-MacLaurin 

series, we obtain 

S = PE-Y.^ ln[l - exp(-/3n^)] = pE + ^ + -\n f3 - - ln(27r) + 0{(3) , (10) 
where 

c(.) = r(i + l)c(i + i/.). (11) 

In the leading order, for determining the stationary point, we ignore the ln/3 term in the 
derivatives of S and keeping only the dominant term we obtain 



Therefore the saddle-point is given by 
The notation may be simplified by setting 



^o^f^l . (13) 



(14) 



\ s J 

s 

so that (3q = Kg E 1+^ . Substituting this value in the saddle point expression for the density 
of states in Eq. © 



pUE) = -^^iTiyJ^E-^ exp [k,{s + l)ETh] . (15) 
(27r)~V* + ^ 

The RHS of the above equation is identical to that given for p'^{n) in P], the number of ways 
of expressing as a sum of integers with s*^ powers, if we replace E by the integer n. For 
s = 1, for example, we have 



exp[7ry 3 
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which is simply the number of partitions of an integer E in terms of other integers. For 
example 5=5, 1+4, 2+3, 1+1+3, 1+2+2, 1+1+1+2, and 1+1+1+1+1, so p(5) = 7. Of 
course, the above asymptotic formula is not expected to be accurate for such a small integer, 
but it improves in accuracy for large numbers. 

While the "physicists derivation" of the number partitions has been known for a while 
and indeed has been extensively used in the analysis of number fluctuation in harmonically 
trapped bose gases the derivation for a general power law spectrum given above is 
novel even though the result was derived long ago by Hardy and Ramanujan p using more 
advanced methods. Equally interesting from the point of view of physics is the sensitivity 
of the bosonic density of states on the single-particle spectrum, in contrast to the fermionic 
Bethe-formula. For example, where as in a harmonic well, both fermions and bosons have the 
exponential square- root dependence in energy for the density of states, as given in Eq. (fT^ . 
in a square-well only the fermions obey such a relation when the low temperature expansion 
is used. For the bosonic case, from Eq. (fT5|) . the density of states is 



This is the same as the asymptotic formula derived by Hardy and Ramanujan for the parti- 
tion of E into squares, for example 5 = 1^ + 2^, 1^ + 1^ + 1^ + 1^ + 1^. It is to be noted that 
in making the identification of p^ln) with p1^{E), E = n is to he identified as the excitation 
energy of the quantum system with a fictitious ground state at zero energy added to the 
square well. 

In Fig. (P) we show a comparison between the exact (computed) p{n) (continuous line), 
and Pooi^) (dashed line), as given by Eq. ()l(ij) . We note that the Hardy- Ramanujan formula 
works well even for small n. Similarly, in Fig. 0, the computed p'^{n) is compared with 
p^(n), as given by Eq. (fT7|) . It will be noted from Fig. (2) that the computed p'^{n) has 
step-like discontinuities, unlike the smooth behavior of p^(n), specially for small n. We 
should remind the reader that these results are not new, and the corrections to the leading 
order Hardy- Ramanujan formula are also known in the number theory literature. We shall, 
however, obtain some new results using our method for distinct partitions (P{n) in the next 
section. 

Before we conclude this section, we note that keeping terms of order /3 in the saddle point 
expansion of S merely shifts the energy E by the coefficient of the term proportional to /3. 




(17) 
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For the s = 1 case in Eq. ()lUp. there is indeed a term hke —^P, leading to the replacement 
of E hj {E — ^) in Eq. (fT^ . The resulting asymptotic expression for the density is the first 
term of the exact convergent series for partitions obtained by Rademacher jl^ . Interestingly, 
for s = 2, there is no term of order (3 in the Euler-Maclaurin expansion. A similar situation 
prevails for distinct partitions as will be shown in the next section. 



ASYMPTOTIC DENSITY OF STATES WITH DISTINCT PARTITIONS 

We now modify the method to obtain distinct partitions of an integer n into s^^ powers, 
to be denoted by d^{n). For example, for s = 1, n = 5, the number of distinct integer 
partitions are 5, 2+3, and 1+4, so d{b) = 3. For distinct partitions, the first guess would 
be to use the fermionic partition function instead of the bosonic one of the previous section 
since distinctiveness of the parts is immediately ensured by the Pauli principle. However, 
there is a problem here which we illustrate using the s = 1 spectrum. For this case the 
fermionic partition function of non-interacting particles is given by (setting x = exp(— /?) as 
before), 

TV 1 oo 

ZNix) = x^'l^ n 7VZ-^. = E ^)^" , (18) 

m=l U ^ I n=0 

which is the same as the bosonic partition function in a harmonic potential, except for the 
prefactor which is related to the ground state energy of particles in the trap. Obviously, 
the f2(A^, ra) is the same for both fermions and bosons even though dj^in) is different from 
•pjsjin). This is because the quantum mechanical ground state of fermions consists of occupied 
levels up to the fermi energy, unlike the bosons which all occupy the lowest energy state. 
Thus, for the fermions at any excitation energy, one should consider the distribution among 
particles as well as holes, each of which is separately distinct and obey the Pauli 

principle. As we show below, the particle distribution at a given excitation energy measured 
from the Fermi energy identically reproduces (the unrestricted) but distinct partitions of an 
integer ra, when n is identified with the excitation energy. 

The relevant "partition" function for the mf spectrum is given by, 

oo 

lnZoo(/?) = E + exp(-/?m^)] , (19) 

m=l 

and the entropy 5'(/3) is obtained as usual by adding (iE to the above expression. Notice 
that this resembles the entropy of an N-fermion system, but with the chemical potential 
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fi = 0. In the normal N-fermion system at any given excitation energy the number of macro 
states available depends on the distribution of both particles above the Fermi energy and 
holes below the Fermi energy in the ground states. By setting /i = we are ignoring the hole 
distribution but only taking into account the states associated with the particle distribution. 
Because of Pauli principle implied in the above form for the entropy, only distinct partition 
of energy E is allowed. Again, using the variable x = exp(— /3) in Eq. (fT^. Zoq{x) above 
is seen to be the genera ting function for distinct partitions d'^{n) of an integer n into s*'* 
powers of other integers 12j |. 

Once this point is noted, the rest of the calculation proceeds as in the case of bosons and 
we obtain the following expression using the Euler-MacLaurin series 

S{f3)=(3E + ^-hn{2) + 0{l3), (20) 



where 



D{s) = ni + -Ul + l/s), (21) 
s 



where 7]{s) = ji — denotes the alternating zeta function. Note there is no log(/3) 
term in Eq. ()2()j) . The saddle point (3o is obtained by setting S'{f3o) = as before. Defining 

A, = {Dis)/sf^'+''> (22) 

and using Eq. (jH)), we obtain 



exp 



'l + s)X,E—s 



PUF)iE) = ^sX, 1 , (23) 



2V7r(l + s)E— 

where the subscript (E) in p is the remind the reader that Fermi statistics has been used 
(with yU = 0). Once again for s=l we recover the well known asymptotic formula for the 
unrestricted but distinct partitions d{n) of an integer ||171], namely 



_ exp[7ry-|] 

P^(F)iE) = 4 ^ 31/4^3/4 ' (24) 

where, as usual, E should be read as n. Similarly the asymptotic expression for d^{n) is given 
by in Eq. (j^Hj) . We have not found this general expression in the literature, though we believe 
it must be known. As we shall presently see, however, a remarkable finding is made when 
this asymptotic formula is compared with the exact computation of d^{n) for s = 2. First, 
however, in Fig. (jS)), we show a comparison of the asymptotic density Poo(F) ^^^1 the exact 



distinct partitions d{n) of integer n for s = 1. As in the case of bosonic partitions p{n), the 
asymptotic formula for d{n) works reasonably, except for n < 10. But the really interesting 
result is shown in Fig. (@]) where we compare Eq. (j23|) for s = 2 with exact computations of 
(P{n). The asymptotic density of states follows the average of the exact (P{n) closely, but 
there are pronounced beat-like structure superposed on this smooth curve. This has come 
about because we have joined the computed points of d?{n) for discrete n's by zig-zag lines. 
Note that compared to d{n), the magnitude of ^^{n) is very small, and this is one reason 
that the fluctuations in d'^{n) look so prominent. We have checked numerically, however, 
that the ratio of the amplitude of the oscillations to its smooth average value decreases from 
about 1.5 to 0.2 as n is increased to 1000. This means that for n — > oo, the smooth part 
will eventually mask the fluctuations. 

Although we cannot analytically reproduce these fluctuations in the many-particle density 
of states (or equivalently in d'^{n)), we can show from the quantum point of view that the 
smooth part 'p'^oo{F) arises strictly from the smooth part of the single-particle density of 
states. To make this point, let us derive the single-particle density of states, g{e), for the 

spectrum. We begin with the knowledge of the exact single-particle spectrum, and write 
the canonical partition function: 

oo 

Zi(/?) = ^exp(-/3n^) . (25) 

n=l 

To express this in a tractable form for Laplace-inverting, we use the (exact) Poisson sum 
formula 

oo oo 

E F{n) = Y: Hq) , (26) 

n=— oo g=— oo 

where 

/oo 
dn F{n) exp{2Tciqn) . (27) 
'OO 



Taking F{n) = exp(— /3?t,^) then gives J-'{q) = y vr/ (/?) exp(— vr^g^/ (/?)). Using this result, we 
obtain 

1 / \ 1/2 / \ 1/2 OO 

^■(«-2(i) -2 + (i) i:exp(-.V/(«) . (28) 

On Laplace-inverting term by term, we obtain the exact result for the single-particle density 
of states : 

11 oo 

9{e) = ^ - 2^(^) + 7^ (2''^^) ' (2^) 

= g{e) + 5g{e), (30) 
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where ^(e) is the "smooth" part consisting of the first two terms on the RHS of Eq. ()29p. 
and Sg{e) denotes the remaining oscillating terms. We can now evaluate Eq. (|T^ for s = 2 
using the above ^'(e) : 

In Zoo = / ^(e) ln[l + exp(-/3e)] de . (31) 
Jo 

Evaluating the integrals, and adding i3E to it, we get the entropy 

We note that the first two terms on the RHS of the above equation are the same as obtained 
earlier in Eq. (j2m) using the Euler-MacLaurin expansion. These yielded the smooth many- 
body density of states given by Eq. on using the saddle-point approximation. The term 
with the double sum in Eq. ()32|) . which arise from ^^'(e) in Eq. (jHUj) and Fermi statistics, 
must be the source of the fiuctuations seen in the density of states in Fig. 4 (the same Sg{e), 
when used in the bosonic case, gives a very different contribution to S{(3)). In principle, 
exact Laplace inversion of exp[5'(/5)], where S{P) is given by Eq. (jH^ . should yield the 
fluctuating degeneracies of the quantum states with E, and hence of d'^{n). We have not 
been able, however, to do this Laplace inversion. Since the oscillation in the exact partitions 
d'^{n) resemble a beat-like structure, at least two frequencies must be interfering to give the 
pattern. Further work is needed to unravel this interesting point. 

FINITE SIZE CORRECTIONS, OR RESTRICTED PARTITIONS 

The smooth part of the many-particle density of states was derived in the previous sections 
for a system with N oo, that corresponded to unrestricted partitions. We now apply the 
same method to obtain the asymptotic density of states for systems with flnite size, that is 
when the number of particles is kept flnite and equal to A^. This corresponds to allowing 
the number of parts to be at most A^. Consider, for example, for s = 1, = 4, n = 5. 
Then, in restricted partitioning, the allowed partitions are 5, 4+1, 3+2, 3+1+1, 2+2+1, 
and 2+1+1+1. The partition with 5 parts, 1+1+1+1+1 is not allowed, since the number 
of parts in this case is greater than 4. The above example is for restricted case that includes 
identical parts in a partition. These will be denoted by p%{n) in general, but for s = 1, the 
superscript will be dropped as usual. For the above example with restricted and distinct 
partitions, however, only 5, 1+4, and 2+3 are allowed. We denote such partitioning by d'lf{n) 
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in general. In this section, we restrict to s = 1, and first present the leading order asymptotic 

expression for pAr(n), using our method of calculating pj^{E). This result is already known 

n 

in the literature by the Erdos-Lehner formula 7], but is derived here because we generalize 
it for obtaining the asymptotic expression for dN{n). To the best of our knowledge, this is 
a new result. 

Asymptotic formula for pN{n) 

The A^— boson canonical partition function in this case is exactly known : 

N 

InZj^iP) = -J2 ln[l - exp(-/5m)]. (33) 

m=l 

The canonical entropy is obtained as before by adding (3E to the above equation. Ex- 
panding the above using Euler MacLaurin series, and assuming that is large so that 
X = exp{—[3N) « 1, even though /3 << 1. We then obtain 

SNif^) = ^oo(/3) - exp(-/3iV)[i - 1], (34) 

The stationary point is determined as before by the condition in Eq. (|7j) and for N large it 
is the same as in Eq. (|T!^ . Substituting this in the saddle point expression for the density 
of states in Eq. (jnj we get 

pAE) =p^{E)exp 

The above expression reproduces the well known correction to the unrestricted partitions 
due to the restriction on the number of particles (see Erdos and Lehner 1^) apart from the 
constant term proportional to ^ in the exponent. This constant may, however, be neglected 
when E is large. Using the conditions /5o << 1 and /3qN >> 1, we see that formula (jH^j) is 
valid in the region C(l) << E << C{1)N'^, where C(l) ~ 1.645. In Fig. 5 we compare the 
two differences, [p^{E) — PNin)], and [pj^{E) — PNin)] for = 20 (Fig. 5a), and = 30 
(Fig. 5b). In the above, p^{E) is obtained from Eq. (fTI)|l . P2o(-^) is the Erdos and Lehner 
formula as given by Eq. (jHSj) . and P2o{n) is the exact (computed) restricted partitions. 
Clearly, the former is much larger than the latter, indicating that Eq. gives a better 
approximation to the exact values for restricted partitions. 




12 



Asymptotic formula for dN{n) 



Next we present the finding of an equivalent asymptotic formula to Eq. ()35|) for the 
restricted and distinct partition. This brings us back to the fermionic particle spectrum as 
discussed in section II and Ref. Eq. (fTUI) of section II does not apply here, however, 

since it is applicable only for the unrestricted distinct partition, ie: N ^ oo. What we need 
is the exact canonical partition function for the particle space. From number theory 16. isj]. 



we found a formula for the (exact) number of ways of partitioning an integer n to at most 
N distinct parts: 

dN{n)=Y.Mn-^^), (36) 

where Pi(n) is the (exact) number of partitions of n to at most i parts, which may be 
generated by the partition function given by Eq. (j33|) . Eq. (j36p implies that the partition 
or generating function for the restricted and distinct partition is given by: 

oo oo i 1 

= n(l + ^")- E U 7T-^y (37) 

n=l i=N+l n=l \^ J 

The first term on the right hand side of Eq. ()37|1 is the generating function for the unrestricted 
distinct partition Eq. (fTIUl . and the second term is a sum of the generating functions for the 
restricted non-distinct partition Eq. ()33p with the integer shifted to i{i + l)/2. To find 
an asymptotic formula for the restricted distinct partition dN{n), as usual, we take inverse 
Laplace transform of Eq. : 

dN{n) = V 11(1 + ^") - E V ^^^n n 



^71=1 J i=N+l L n=l 

oo 



X") 



din) - E Pi('^ ~ ^) ' 

i=N+l 

oo 

PooiF){E)- E UE-A) 



i=N+l 



exp[7r^y] 

4 X 3V4E3/4 ^^^^^ 



^ p^(E-A)exp 



/Q{E-A) 




TT 

= PN(F)iE) . (3 

where A = i[i + l)/2, x = exp{—(3) and n is identified with E. Note that since the 
asymptotic expression for the restricted partition 'p^{E) is valid only for C(l) << E « 
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C{1)N'^, Eq. ()38|) is thus valid only in this range. Fig. 6 displays the two differences, 
p^(^p^{E) - dN{n)], and [p^(^)(^) - d^vH] for = 20 (Fig. 6a), and = 30 (Fig. 6b). 
In the above differences, p^(^p-^{E) is obtained from Eq. (f^ . pj^(^p-^{E) from Eq. (jSHj), and 
dN{n) is the exact (computed) restricted distinct partitions. Again, similar to the non- 
distinct case (Fig. 5), the N-correction asymptotic formula gives a better approximation to 
the exact finite N partition than the infinite distinct one. 



DISCUSSION 



This work emphasizes the connection between the many-body quantum density of states 
in a power-law spectrum with the number theoretic partitions p'^{n), and the distinct parti- 
tions d^{n). This was already well-known to the physics community for p{n). Most of the 
results derived in this paper are known in the mathematical literature, with the possible ex- 
ception of the the asymptotic formula for d^{n) (Eq. ()23|l ). and the generalized formula ()38|l 
for restricted distinct partitions. Perhaps the most interesting result in the paper is shown in 
Fig. 4, where the fluctuations in d'^{n) are shown. These oscillations in this number-theoretic 
quantity are linked, from the quantum mechanical point of view, to the oscillating part of 
the density of states in a square-well potential, and the Pauli principle. However, we are 
not able to completely demonstrate this point to our satisfaction because of the difficulty of 
Laplace inversion of exponentiated quantities. This is left for the future. 
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FIG. 1: Comparison of the exact p{n) (solid line) and the asymptotic p^{E) (dashed line), obtained 
from Eq. H16|) for s = 1. 

FIG. 2: Comparison of the exact p^in) (solid line) and the asymptotic p'^{E) (dashed line), 
obtained from Eq. H17|l for s = 2. 

FIG. 3: Comparison of the exact d{n) (solid line) and the asymptotic Poo(F)(^) (dashed line), 
obtained from Eq. for s = 1 and distinct partitions. 
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FIG. 4: Comparison of the exact cf{n) (solid line) and the asymptotic 'p'^^i^p^{E) (dashed line), 
obtained from Eq. H23j) for s = 2 and distinct partitions. Note that the y-axis is no longer in log 
scale. 

FIG. 5: (a) Comparison of [Poo(-^) ~J'2o("-)] (dotted line) and \p2o[E) — P2o('^)] (solid line) for 
N = 20, where 'p^{E) is obtained from Eq. (jlbf) . P2o(^) the Erdos and Lehner formula as given 
by Eq. (jHSJ, and P2o{n) is the exact (computed) restricted partitions, (b) Same for N = 30. 



FIG. 6: (a) Comparison of Poo(_F)(^) ~ d2Q{n) (dotted line) and P2o{F)iE) — d2o{n) (solid line) 
for = 20, where Poo{F){E) is obtained from Eq. ([21]). P2o(F){E) from Eq. (|38|) . and d2o{n) is the 
exact (computed) restricted distinct partitions, (b) Same for = 30. 
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